in 



(N 



X 



Analysis of data 
on low energy nN ttttN reaction. 
I. Total cross sections. * 

V.V. Vereshagin", S.G. Sherman^, A.N.Manashov", 
U. Bohnert^ M. Dillig^, W. Eyrich^ O. Jakel'^, M. Moosburger^ 

") St. Petersburg State University, St. Petersburg, Russia 



^ '. ^) St. Petersburg Nuclear Physics Insitute, St. Petersburg, Russia 



^) Physikalisches Institut IV der Universitat Erlangen-Niirnberg, 

Erlangen, Germany 



\0 • J. . .. 

cn ■ jinstitut fiir Theoretische Physik III der Universitat Erlangen-Niirnberg, 

^ ; Erlangen, Germany 

February 1, 2008 



Abstract 



This is the first of a series of papers on a consistent model independent analysis of 
I the complete experimental information on the reaction vr — > vrvr at pion momenta 

up to 500 MeV/c. The paper summarizes the theoretical approach and details of 
the computational procedure. The complete database on total cross sections in 5 
ttttN channels is given together with a critical discussion of their model independent 
analysis. 
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1 Introduction 



The ttA^ — i> TinN processes near threshold is currently a subject of intensive investiga- 
tions, both theoretically and experimentally [1-25]. There are two main reasons why this 
process has attracted growing interest: a) the development of experimental techniques 
made it possible to measure in the nearest neighborhood of threshold not only total cross 
sections but also angular distributions; b) reliable data on the parameters of low energy 
TTTT interaction are necessary for a precise estimate of higher order corrections in Chiral 
Perturbation Theory (ChPT), which gives a link between QCD and the world of hadrons 
(see, for example, [^). The reaction in question provides an excellent opportunity to get 



the parameters appearing in the ChPT expansions |2^, 29 . 

In view of the importance of reliable and sufficiently precise results (which has been 
specially stressed in a recent paper ||3^, devoted to the problem of estimating of higher 



order terms of ChPT) it is desirable to perform a complete, model independent analysis 
based on the full set of available experimental information on vrA^ ttttN reaction near 
threshold. The main goal of such analysis would be to establish the phenomenological 
form of the reaction amplitude in the low energy domain and to fix (with the help of data) 
the values of the corresponding parameters appearing in this form. 

It is clear that this rather complex task is preferentially solved in various steps, as 
the systematic experimental and, especially, theoretical investigation of the low energy 
dynamics of the reaction in question is, in fact, just beginning. The problem of the 
extraction of low energy vrvr - interaction parameters with a desirable accuracy looks, in 
such a situation, extremely difficult. At the same time, the successive resolution of this 
problem is important by its own, irrespective of the following application of the results: 
the pion, which is the lightest hadron, participates — really or virtually — in every strong 



interaction process (see, for example, |£8|, ^ 531 )• This is a reason why the complete 



phenomenological analysis of the data on low energy ttN tttcN reaction is really of 
actual interest. 

In the present paper we make a first step in this direction. It combines two points: it 
includes the phenomenological ansatz for the amplitude (the generalization of the standard 
approach used in the isotopic analysis) and the corresponding analysis of the full set of 
experimental data on total cross sections of 5 channels at Piab < 500MeV/c. The goal of 
this paper is to answer the following questions: 

1. Is the modern database on total cross sections internally consistent or are there in 
a statistical sense "doubtful" points? 

2. How strong is the comparative infiuence of the various mechanisms (tt - and A- 
exchanges, "background") on the vrA^ -ktiN reaction amplitude? 

3. What is the capability of data to distinguish between the different contributions? 
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4. Is it possible to extract the values of the vrA^ ttttN amplitude parameters from 
the data on total cross sections only (of particular interest being the values of the 
parameters describing the low energy vrvr-interaction)? 

The new precise data on vrA^ tcttN total cross sections, which became available in 
the course of the last years (see Table 1) support the common expectation to contain 
accurate information on the low energy vrvr interaction near threshold. As pointed out 
above, this information would be of great importance in connection with the problem of 
a quantitative description of the low energy hadron dynamics from QCD. 

To extract the information on tttt - scattering from the data on ttN ttttN reactions, 
one has to perform a refined analysis, because the detailed mechanism of the latter reac- 
tion is not understood in detail and the influence of the 47r-vertex being not necessarily 
dominant. Its contribution could be (if at all) separated among all others with the help 
of an accurate analysis of the full set of experimental data including those on differential 
cross sections. 

At this moment, however, the overwhelming majority of recent results on tttt- 
parameters is based on analyses of data on total cross sections only. Most of the analyses 
having been performed with the help of the Olsson-Turner (OT) model [EH] (for a review 



see [0). This situation is similar to that of the seventies (see, for example, refs. ||2^, pH]). 
That time the most popular method of the corresponding data analysis was the Chew-Low 
extrapolation procedure This method, however, proved to be not sufficiently reliable, 
the main uncertainty originating from the phenomenological features of the vrA^ —>■ ttttN 
reaction: the contribution of the OPE graph turned out to be small in comparison with 
those of the other graphs (see, for example, the recent papers |^8|, |59|, Q). Due to this 
reason the results of the extrapolation of experimental data to the unphysical point r = /i^ 
become extremely sensitive to small variations (errors) of the data used. From the purely 
mathematical point of view this effect is just a manifestation of J. Hadamards general 
theorem concerning the ambiguities of extrapolation problems (see ref. ||33|| ) . 

The same situation repeats now again. A good example is given by the comparison 
of the results on S-wave TTTr-scattering lengths obtained by different groups, where the 
most interesting conclusion has been made by the OMICRON group in ref. @]: a single 
parameter ^ of the OT model is not sufficient for the adequate description of the low 
energy data on total cross sections of three channels {n~n~^n, 7r~7r°p, vr+Tr+ra). It should 
be stressed here that from the purely technical point of view the analysis made by this 
group is one of the most accurate, because it accounts for all specific features of the OT 
model. It looks much as if the origin of the difficulty mentioned above is hidden either in 
these very specific features of the OT model or in the inconsistency of the data. 

The original formulae given in ref. |^ are valid only for the threshold values of the 
amplitudes. Thus, to extract the value of ^ (or, equivalently, S-wave scattering lengths) 
one has to perform the extrapolation of the experimental data on total cross sections to 
the threshold in all 5 channels of the ttN — ^ TrvriV reaction. On this way, however, one 
meets at least two serious obstacles closely connected to each other. First: to reduce 
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errors one must use for the extrapolation the datapoints nearest to threshold. Second: 
one needs to know the form of the extrapolation law. As it is well known the low energy 
measurements exhibit extremely delicate experimental problems by itself, the difficulties 
increase approaching to threshold. Thus, those very points which would be the most suit- 
able for the extrapolation, are the most difficult to measure. It must be mentioned that — 
in addition to purely experimental difficulties — there is the unresolved theoretical prob- 
lem of data interpretation: in the near threshold domain the radiative corrections become 
significant. Furthermore, the form of the extrapolation law is more or less known for the 
case of elastic scattering only (due to general requirements given by quantum mechanics, 
such as partial wave expansion, unitarity, etc.). In the case under consideration — the 
vrA^ — s> vrvrA^ reaction — these limitations are considerably weaker and the extrapolation 
law itself is still unknown. 

All these obstacles result in large systematic errors of the parameters to be determined. 
Besides, needless to mention, the OT model — as every other one — has its own limits 
of applicability. 

Out of this situation we started a general analysis of the phenomenological features of 
the total cross sections of the 5 tiN ttttN channels in the low energy domain {piab < 
500MeV/c). The main goal of this analysis is to get a deeper understanding of the 
phenomenology of the process and to estimate the informative capacity and selfconsistency 
of the modern database. 

This paper is the first one of a series devoted to the systematic model independent 
analysis of all available experimental information on the low energy tcN — > ttttN reaction. 
It contains — along with the analysis of data on total cross sections — the detailed 
description of the theoretical approach and the computational procedure, as both subjects 
are important for the subsequent analysis of distributions also. That is why the text is 
divided into 4 sections. The details of the approach are discussed in Sec. 2, where we 
present also the complete database on total cross sections. Sec. 3 contains the results of 
the theoretical analysis of the experimental situation and those of the data fitting . The 
discussion of the results obtained and the concluding statements are given in Sec. 4. 

2 Theoretical approach 

2.1 Decomposition of the amplitude 

The main ingredients of our approach are given in refs. 0, ^, |19|. Here, we just want to 
recall some important points and give an extension of the formalism. 

Our ansatz is model independent in the sense that we do not make a priori suggestions 
on the form of interactions (other than the well known symmetry requirements), coupling 
constants (except the nNN coupling) or on the off-shell behaviour of the different ver- 
tices. It is based on the exact symmetry principles (Lorentz-invariance, isotopic and 
crossing symmetry, C-, P- and T-invariance and Bose-statistics for pions) and some well 
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established features of the low energy phenomenology, the basic idea being precisely what 
is used for the solution of a so called "incorrect" tasks (frequently also called "ill-posed" 
or "ill-conditioned"; for a detailed review of this subject see ref. |Q and refs. quoted 
therein). Our procedure can be schematically formulated in the following way: 
1. Write down the amplitude M of the process under consideration in the form: 

M = S + B, (1) 

where S stands for the "special" contributions with a strong dependence on kinematical 
variables, which must be separated and where B includes all other contributions (in 
accordance with the terminology of ref. we call B the "background"). 



2. Take into account all available information on the form of the functions 5* and B. 
In our case, this step implies the selection of convenient parameterizations for the above 
mentioned functions in terms of the kinematical variables suitable for our problem. Thus, 
one obtains the amplitude M of the reaction in question as a function of some set of 
free parameters. From the decomposition of the invariant amplitude M — as motivated 
by the underlying dynamics — into S*, with a sensitive dependence on variables (as will 
be discussed below, S contains the OPE graph and the A-dominated amplitudes), and 
a smooth background S, it is clear that the parameterization of S and B differs in its 
functional structure. 

3. Fit experimental data with the amplitude M, obtained in step 2, to determine the 
influence of the various elements of M on the data and to find, if possible, the values of 
the parameters or at least some combinations of them. 

It is clear, that the more detailed information on the form of the functions 5* and B 
is involved, the stronger is the model dependence of the results obtained. Hence it is 
suitable here to fix the meaning of the term "model independent approach". To do this, 
we give here in a shortened form the necessary formulae (detailed explanations are given 
in ref. 0). 

Ti,{h) + N^P) ^ vr,(fc2) + T^c{h) + iV;(g) (2) 

(a, b,c = 1, 2, 3 and /, Z = 1, 2 are isotopic indices; /x, z/ = 1, 2 are the nucleon polariza- 
tions), can be written in the "canonical" form: 

M;f^^ = (r°)ji(5^^A'^^ + {T^)fi6''''B>"' + {T^)fi6''''C^"' + te^^'6fiD>"' , (3) 

where {X = A, B, C, D) 

= M'^(g)Xu^(p), (4) 

X = {^,.J + v;.fc + v;-^ + ^T,[A;,l]}(275), (5) 

k = -ki + ek2 + efcs, k = -ki + ek2 + efcs, (6) 

e = exp{i—), e = e* , (7) 



The amplitude M'i\'^ of the reaction 
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and K = 'j^K^^ with K = k,k. The 16 invariant form factors Sx,Vx,Vx,Trc (with X = 
A, B, C, D) are functions of five independent scalar variables: 

r = (p-g)2, iy={p + q)-k, e = {p-q)-k, 

u = {p + q) ■ k, 6 = {p — q) -k. (8) 

The general properties of these functions with respect to the requirements of Bose, 
crossing and C-symmetry are described in |[T^ and references therein. 



So, up to this point, no suggestions have been made on the form of the amplitude M 
— only the restrictions due to the exact symmetries of strong interactions are being taken 
into account. 

To proceed further, however, one has to include some additional physical information. 
First of all, it is necessary to choose the energy domain, which is the most suitable for 
the solution of our task — the determination of vrvr-scattering parameters. 

In a previous paper we have discussed only general features of our method. The 
particular form of the parameterization considered there is suitable only for a sufficiently 
narrow interval of the incident beam momentum ki. In such a case the possible influence 
of resonances (which were implied as being far from the boundaries of the physical phase 
space) can be described by slowly varying real functions. 

Thus, in a narrow momentum range, physical arguments suggest a representation of 
M as a sum of two terms: 1. the contribution of the OPE graph (denoted hereafter as 
n); 2. the total contribution of all graphs besides OPE ("background" denoted as B). 
The crucial point is to note that, inside of this narrow kinematical interval, the first term 
(OPE) is characterized by a strong r dependence (due to the pole which is placed close 
to the boundary of the physical phase space), while the second one ("background") has 
no such strong dependence on the variables because of the absence of the corresponding 
physical mechanisms. Reversely, this means that the background contribution can be 
approximated by any smooth function of the variables from eq. (H). The choice of the 
approximating functions is certainly not unique (this weakens the meaning of the term 
"model independent approach"). 

In this note — compared to [|1^] — our goal is more ambitious; we would like to extend 



our analysis to the full kinematical region from threshold (or some MeV above it, to avoid 
problems caused by the rescattering of final particles; see, for example, [^) up to ki iab ~ 
0.5 GeV/c. In this domain the vrA^ tittN reaction amplitude is strongly infiuenced by 
the A-isobar creation in the initial {ki^iab ~ 0.3 GeV/c) and final {ki^iab ~ 0.5 GeV/c) ttN 
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In this case the contribution of the graphs with A isobar exchange (simply called A 
below) shows a sensitive dependence on the scattering energy due to its resonant structure 
in the initial and final vrA^ invariant masses. Therefore, it has to be separated from the 
background amplitude. This contribution cannot be approximated by linear functions in 
the kinematical variables with the desirable accuracy. Furthermore, the behavior of the 
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amplitude near the resonance position is strongly restricted by unitarity. Consequently, 
our amplitude has to be extended: we must add the A-isobar contribution explicitly, for 
which, in particular, the presence of an imaginary part is necessary. 

Keeping in mind all these notes above, we arrive at the conclusion that the appropriate 
way to construct a suitable approximation for the amplitude M — in the momentum range 
investigated — is the decomposition: 

M = B + U + A (9) 

(i.e. = n + A, eq. (|l]), see also fig. 1), with the following amplitudes: 

- B is the total contribution of all graphs besides those with A poles and OPE; It can 
be approximated by linear polynomials with unknown coefficients. This procedure 
gives 11 free parameters: Ai - An, (see Sec. 2. 2.1 below). 

- n is the main ingredient of the picture — the OPE graph. It contains four free 
parameters: Qq - g-^, (see Sec. 2. 2. 2 below). 

- A is the A-isobar contribution which can be approximated with functions of the 
Breit-Wigner type (with the terms zMaPa in the denominator taken into account). 
The simplest way to write down the unknown polynomials in the corresponding 
numerators is to use effective Lagrangians, this method guaranteing the validity of all 
necessary symmetry conditions (isotopic and discrete symmetries, etc.). According 
to Sec. 2.2.3 (see below) this introduces 6 extra free parameters. 

Taking everything into account, one obtains the amplitude M of the reaction (1) in 
the form given by the formulae (^|-|^), with the functions Sx, Vx, Vx, Tx {X = A, B, C, D) 
fixed by 21 free parameters to be determined with the help of experimental data at 
ki,iab < 500 MeV/c. 

We have to remember that only 10 of these 21 parameters (namely, those, describing 
the OPE and the A contributions) can be interpreted as coupling constants of the cor- 
responding vertices, while the 11 remaining in the background amplitude have no clear 
relation to underlying physical parameters (e.g. the iV*(1440) position and width). 

Some of 10 coupling constants mentioned above could be fixed (or, at least, con- 
strained) by additinal theoretical and experimental information. For example, the vrA^A 
coupling could be fixed by the known value of the isobar decay widths (or, equivalently, by 
the additive quark model) while the AvrA one - by the requirements of SU{Q) symmetry. 
We, howeveer, prefer - for the sake of fiexibility of our procedure - to keep them as free 
parameters. This allows us to avoid the introducing of unnecessary model dependence 
into our approach. 

Of course, it is not to be expected that the restricted data set on total cross sections 
only is able to fix all 21 parameters in quantitative way. To solve this task one has to 
analyze the full set of distributions at several values of incident momentum. Our present 



6 



goal is rather to answer the questions hsted in Sec.l and to give insight on degree of model 
dependence of the previously known results. 

In closing this section we remark that in order to explore the significance of the various 
pieces of M, both fits with M = B,U,A only and with M = U + B,A + B,U + A were 
performed in comparison to the fit with the full amplitude M = B + H + A. 



2.2 The amplitude 

The explicit inclusion of the A contribution is the main extension of the formulae given 
in ref. Therefore we give in the following a detailed description of the formal- 

ism; for completeness, we supplement this presentation by the baswic ingredients for the 
parametrization of the OPE and the background terms from ref. . 



2.2.1 The background contribution 



As shown in the aggregate contribution of the backgroud mechanisms can be described 
as follows: 



Sa = Ai + A2T + A3{9 + 9) 

Va = A4 + A^T + AfiO + ArO (10) 

Va = Ai + A5T + AqO + A^e 

Ta = lAsii^-u) 



Sb = Ai + A2T + As{e9 + e9) 

Vb = A^ + eA^T + AqO + eAjO (11) 

Vb = eA4 + eA^T + AqO + eAjO 

Tb = tAs{eu - eu) 



Sc = Ai + A2T + Asiee + eO) 

Vc = eA^ + eA^T + A^e + eArO (12) 

Vc = eAi + eA^T + AqO + eA-rO 

Tc = iAs{eu - eij) 



Sn = 
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Vd = lA^u (13) 
Td = Alo + AuT 

Here the parameters Ai -^ An are real numbers and the kinematical variables 
T,nu,h',9,9 are defined by eqs. (||) above. 

2.2.2 The OPE graph contribution 



According to ref. [T^ the OPE graph contributes to 5*^, 5*^, Sc functions only, the corre- 



sponding expressions having the following form. 



Sa = 




2g 


T 


- 


Sa = 






T 


- 


Sa = 






T 


- 



-(Go + Gi{e + e) + G2{ee) + G^{e^ + e^)) 

-(Go + Gi(e0 + tO) + G^m + G-i{te^ + W^)) (14) 
-(Go + Gi{ee + W) + G2(#) + G^ieO^ + e^')) 



Here Gq, . . . , G3 stand for free real parameters describing the vrvr interaction and the 
constant 9 corresponds to the ttNN vertex. It can be expressed via the familiar pseu- 
dovector coupling constant / {p/An ~ 0.08) as follows 

2m 
9 = / 

2.2.3 The A contribution 

In accordance with sec. 2.1 the propagator of the A-isobar involves poles of the Breit- 
Wigner type: 

ul^-Ml + tM^TA' ^^^^ 
Here, uj^^n stands for the invariant mass of the ttN pair. There are 6 such pairs: 

u^li (^ = 1,2,3): {p + hf, {p-k2f, {p-h)\ (16) 
(z = 1,2,3): {q-h)\ {q + h)\ {q + h)\ (17) 

The numerators of the corresponding terms must have the correct spin-isospin struc- 
ture (Ja = /a = 3/2); they have to fulfill also the correct properties with respect to Bose 
and crossing-symmetry requirements. The simplest way to coordinate all these conditions 
is to use the method of effective Lagrangians, the corresponding coupling constants being 
regarded as free parameters. 
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Though we expect a pronounced energy dependence (A-pole) only for the direct graph 
of fig. 2d (the additional nucleon - and A -poles appearing in graphs of figs. 2a-2d cannot 
give rise to the strong dependence of the kinematical variables because they are located 
too far from the physical space) we include in our analysis all graphs of figs. 2a-2d (along 
with crossing ones) for the sake of consistency. 

To calculate the contributions of graphs of figs. 2a-2d we use the following effective 
Lagrangians: 



C^NN = gnNNNr^^l^l^Nd^TT'', (18) 

jC.na = g^NA^lNdV + h.c. (19) 

= g^AA^Itl^A^'^^'^'n^ (20) 

= N{hF^'' + ?>hF^'')ii,^ld^'K\' + h.c. 

+ iV((7oFo"''' + 3^7ii^f^)7.75A^5'^7rUvr^ + /i.c., (21) 



where 



pabc ^ i^abc ^ ^ab^c _ ^ac^b ^22) 

2 
3 

are the invariant isotopic forms guaranteing the correct isospin (1=3/2) of the isobar field 

r^Aa = (23) 

(r° are the standard isospin matrices: [T",r^] = ie'^^'^t'^). 

It should be stressed here that the choice (pO|) of the vrAA coupling (non-derivative 
one) has nothing to do with the problem of chiral symmetry breaking, because (due to 
the equation of motion) the chiral invariant form of the ttAA coupling 



is equivalent to that of the eq. ( pil| ) on the mass shell of the isobar. The off-shell 
contributions renormalize the constants of the background and the ttttNA vertex. Since 
the parameters of the background as well as those of the A-vertices are considered to be 
free, the difference between "chiral" and "non-chiral" types of vrAA coupling becomes 
unimportant; our choice of the vrAA vertex is dictated by simplicity. The same is true 
for the ttNN vertex also. We use the familiar pseudovector type ([T8|). 

This note also explains the absence of the terms with the derivative of the nucleon 
field in (^TJ) — such a term could again only renormalize the parameters already taken 
into account. 
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The effective Lagrangians (p!8|-pTD describe the A contribution to the vriV — > -mxN 
amphtude (in the energy region under consideration) in the most general way. They 
contain 6 real free parameters {gnNA, Q-kAA, fo, fi, Qo, di), one of them {gnNA) can be fixed 
in accordance with the known A — >■ ttN decay width. We leave it free due to the reasons 
explained in sect. 2.1. 

To complete the description of the formalism we give the form of the A propagator 
used in our calculations: 



^-(^) = (24) 



Pm7^ + Ma 

P^-MI + iMaTa 

2 1 



M and F being the isobar mass and A ttN decay widths respectively. 



2.3 The fitting procedure 



To fit experimental data one has at first to calculate the modulo squared | M | (averaged 
over initial and summed over final nucleon polarizations) of the theoretical amplitude for 
every channel of the reaction (1) at any fixed values of the free parameters Ai {i = 1, 21). 
The second step consists of the integration of | M p over the appropriate part of the phase 
space, the total energy E being fixed. This gives the theoretical expression: 

^'^'^^ 'E=const= V A;;,gW(E,xi,...,x,)A„, (l<r<4). (26) 



w 

m,n=l 



dxi...dx. 

The distribution V to be fitted is given by a convolution: 

V^'\E,Xi,...,Xr) = ^ '^{E,Xi,...,Xr)*<^{E,Xi,...,Xr), (27) 

OXi...OXr 

where $ is a known apparatus function specifying the type of distribution, experimental 
conditions, etc.. Expression ( p6|) is a bilinear form in the free parameters Ai with the 
coefficient matrix being a function of the set of variables {E, xi, ...,Xr) describing the 
distribution under consideration. 

To perform the fitting procedure, it is necessary to organize the experimental distri- 
butions in the form of a finite array of numbers. To do this one has to divide the phase 
space into a set of N cells Cj with "central points" Pi{E, x\, x*), with i = 1, ...,N. Each 
point Pi must be supplied with the corresponding experimental number Wi (integral over 
the volume of Cj) and, of course, with the suitable error bars. In other words, one obtains 
a r-dimensional histogram representing the experimental data in question. 
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Precisely the same procedure must be performed for the theoretical distribution on 
the right hand side of (pBf ). This means that for each central point Pi of the distribution 
one must prepare the corresponding 21 x 21 matrix Q^^'^{E, x\, ), N matrices in 
total. When several different experimental distributions are being treated simultaneously, 
the necessary number of correlation matrices increases proportionally. Finally, to treat 
data at different energies, one has to prepare the above mentioned full set of correlation 
matrices for each energy value. Keeping this in mind, one can easily imagine the size of 
computer memory which is necessary to perform the data treatment. (This also explains 
why our program could not be run earlier.) 

After all preliminary steps have been done, the fitting procedure is performed with 
the help of the standard program FUMILI . 



2.4 Integration over the phase space 

The question of the integration over the phase space is extremely delicate. The crucial 
point is that near the threshold the phase space factor strongly depends on the values of 
particle masses (this sensitivity was pointed out recently by R.E. Cutkosky — see citation 
in ref. |2T|). The dependence of the amplitude modulo squared is weaker (corrections are 



proportional to mass differences), but, nevertheless, it must be taken into account when 
one uses the standard version of the Monte Carlo program FOWL (or any other 
routine) for calculation of the integrals. 

Here, we specify the remark by R.E. Cutcosky and show that the amplitude modulo 
squared contains two different parts, one of which requires special caution with respect 
to the values of masses used in the course of computation. To explain this statement we 
must return back to expressions (^) and (^. According to these formulae, the amplitude 
modulo squared can be written in the form: 

I E I x^"" \'= E ^;m,,x„ (28) 

^J.,u=l i,j=i 
where X = [S, V, V, T) and the matrix M is hermitian: 

Mi^ = ^Tr{ri{q + mf)r^{p-mi)}; 

T = {I,kX'^[kM, t = {lXk,l[kk]}, (29) 

with f = 7^r'^ for r = p,q,k,k; rrii and rrif are the masses of incoming and outgoing 
nucleons, respectively. 

The matrix M degenerates at the boundary of the phase space and the quadratic form 
(p8|) loses the property of positive definiteness outside of the phase space. The bounds 
of the phase space domains determined at the physical values of masses rrii rrif and 
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at the average ones rrii = mf = rfi do not coincide. This means that one must use the 
physical masses along with momenta to compute both the phase space factor and the 
matrix elements Mij. Otherwise, incorrect and even negative contributions to the left 
hand side of ( ^Sf ) are unavoidable. The point is that errors must be compared not with 
the particle masses, but with the distances from the boundaries of the real physical phase 
space, the latter being small if the energy is low. 

It is easy to understand that all above arguments are valid for the treating the mass 
differences of the pions also. 

Thus, we conclude that one must use the true mass values both for the phase space 
factor and for the matrix elements ( p9D to avoid problems with the kinematics near thresh- 
old. One can take the averaged values m for the nucleon and fl for the pion masses for 
the computation of the amplitude vector X in (|28|); in this case the errors will be of the 
order of the isospin symmetry breaking. 



2.5 Description of the database 

In our analysis we use the full set of published (up to May 1994) experimental results on 
total cross sections of 5 channels (namely, 7r~7r+n, ir'n^p, 7r°7r°n, 7r"'"7r°p, Tr+vr+n) of the 
7tN — >■ irirN reaction at Piab < 500MeV/c. These data (105 points in total) along with 
the corresponding references are listed in the Table 1. For convenience the same data are 
presented graphically on Figs. 3a - 3e in the form of the r-th (r=l,...,5) channel quasi- 
amplitude Mr dependence on the incident beam momentum piab- The quasi-amplitude is 
defined as follows 

= {im^pi^y/^ ^ {^}V2 (30) 

(see Table 1), where dTs is the standard element of the 3-particle phase space volume. 
This method is much more convenient than the one conventionally used {atot as a function 
of Piab), because the extremely rapid energy dependence of the numerator caused by 
the phase space volume is cancelled by the denominator. The resulting smooth energy 
dependence of the quasi amplitude Mj provides a fairely direct insight in the quality of 
the various data points. As an immediate consequence the point no. 12, being obviously 
inconsistent, was excluded from the data set analysis from the very beginning. 



2.6 Strategy of the fitting 

To answer the questions raised in sec. 1, we choose the following procedure: first, we 
fit the total cross section data in order to see what we can learn from this information 
alone. To test the sensitivity of the data to the OPE or A contribution, we begin to fit 
with the background amplitude B only (11 parameters). For comparison with models it 
is interesting also to fit these data using the OPE amplitude only (4 parameters) and, 
for completeness and to test the sensitivity, using the A amplitude only. At this point 
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it has to be stressed, that the amphtude A alone is not the full amplitude which can 
be calculated from the diagrams, since its smooth part corresponding to the off-shell A 
interaction is included into B. In contrast, the amplitude 11 itself has a clear physical 
meaning because it corresponds to the simplest variant of the OPE model (see ref. [|T9l). 

More physical, however, are the results of a fit with a combination of two of these 
amplitudes, like B + A or B + U. At the end, the fit including the full amplitude with 
21 parameters is performed. The full amplitude, according to sec. 2.1, should be general 
enough to describe every reasonable collection of data points in the momentum region 
in question. Any set of data points tolerating a reasonable fit can therefore be considered 
as being selfconsistent. 

The quality of a fit is considered to be good if the value ^ 1, where 

1 ^ 

f^^T^xl (31) 

i=l 

with A^ being the numer of datapoints, whereas M being the number of degrees of freedom 
(A^ minus number of fit paramaters). When there are some data points with individual 
contributions of xf more than 5 (which is of course our own, arbitrary upper limit), these 
points need a closer look and have been excluded for the next step of the fit. 

After the exclusion of some data points the fit with the full amplitude and some fits 
(leading to a large before) for a reduced set of parameters are repeated to see the 
influence of these data points on the results. 

3 Analysis of data on total cross sections 

As mentioned above our analysis of the vrA^ titiN reactions is based mainly on the 
results from ref. supplemented by the inclusion of the A-isobar. Even in this fairly 
restricted parametrization the model - as reflected by its dependence on 21 parameters 
- is highly complex; a finding, which is stressed by the result of the fitting procedure 
from section 4, that none of the various ingredients from eq. (9) of the model strongly 
dominates the complex dynamics at low energies. Due to this subtle interplay we expect 
that the various vrA^ ttttN isospin channels exhibit quite different features, such as in 
their extrapolation to the corresponding threshold or in their angular distributions. To 
gain deeper insight into our subsequent results and to facilitate their interpretation, we 
study in appendix A "quasi" vrA^ elastic scattering (i.e. equal nucleon and pion mass) as 
a much simpler system with related dynamics. This consideration will help to understand 
better the logic of the subsequent analysis of the vrA^ vrvrA^ reaction. 

Our analysis of the vrA^ ttN reaction at low energies exhibits the following conclu- 
sions: 

1. Dealing with reactions of particles with spin one has to take account of the exis- 
tence of several independent scalar amplitudes. In such a case the term" simplest 
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dynamics" loses its transparency and must be specially defined. 

2. The spin structure of the amplitude manifests itself even in the case of unpolarized 
particles: it causes the " minimal" (unremovable) dependence of the distributions on 
kinematical variables (due to the presence of matrix G — see eq. ( |A.7| )). It can be 
shown that it is impossible to construct a physically correct amplitude which would 
result in the precise "phase space-like" behavior of differential cross-sections in all 
channels ( |A.5D . 

3. The explicit form of the "minimal" kinematical dependence mentioned above is 
determined by the interplay of the values of parameters describing the independent 
scalar amplitudes Ai. It may be rather complicated even in the "simplest" cases 
when all Ai are constants. In what follows we call this phenomenon the "spin 
correlation" . 

4. The isospin conservation law connects, to some extent, the spin correlations in dif- 
ferent channels. Then, crossing symmetry and Bose statistics give the additional 
limitations on the admissible form of the latter correlations. Altogether, these sym- 
metries result in the existence of spin-isospin - Bose - crossing correlations con- 
necting the kinematical dependence of the amplitudes of different channels. It is 
easy to understand that the more identical particles are involved in the reaction the 
stronger are these correlations. 

5. The low energy behavior (the form of the extrapolation law) of the total cross section 
in a given channel depends on the relative magnitudes and signs of the amplitude 
parameters; it may be different in different channels (see Appendix, Ex.III). To get 
the form of the extrapolation law one has to analyze the data on distributions at 
several values of incident momentum (see Appendix, Ex.I-III). Example IV shows 
that the "natural" terms (that of order k^) may be absent in low energy expansions 
of total cross sections of all channels. 

Keeping in mind these basic findings from the the "quasi" vrA^ elastic scattering we 
can start now with the analysis of the modern database on total cross sections of the more 
complicated TriV irirN reaction in the low energy domain P^. < 500MeV/c 



3.1 General theoretical consideration. 



First of all, let us summarize some formulae concerning the vrA^ 

m 



we use the notations of ref. 



ttttN reaction (hereafter 
The five experimentally accessible channels of this 



reaction along with the corresponding amplitudes and statistical factors fr are listed 
in Table 2 (in comparison with ref. |]19| the numbers of two last channels are changed). 
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In accordance with the eq. ( pQ]) the total cross section of the r-th channel can be 
written in the form 

ar = N j duj ■ fr ■ Xr^MXr (32) 

Here, is the normalization factor (see eq.(57) of ref. [|19|), X^ — the amplitude vector 
of r-th channel, i.e. 

Xj. = (^Sj-, Vri '^ri Ty) (33) 

and the hermitian 4x4 matrix M (which is the same for all channels, if one neglects the 
isospin symmetry breaking) has the elements given by eq. (pOj) . 

With the help of eq. (p^ it is easy to analyze the correspondence of some statements 
concerning the dynamics of the reaction in question with the experimental data on total 
cross sections. 

I. The OPE - dominance hypothesis. 

According to popular wisdom the main contribution to the low energy -kN ttttN 
amplitude originates from the OPE-graph (11; fig.l). It is easy to check that in this case 
(we assume the leading order of 47r-interaction and the arbitrary form of the nucleon form 
factor) the amplitude of the reaction (|^) takes the form: 

Sa = Sb = Sc = S{t), Sd = 

Vx = Vx = Tx = 0, {X = A,B,C,D;eq.{3)) (34) 

Here, S{t) depends solely on r (see eq. |]). Using the expressions from Table 2 and 
eq. (|32D , one obtains the following relations among the total cross sections of five basic 
channels: 

1 1 

-cr(7r^7r^r?,) = a(n^n^p) = a(n'^n^n) = a(n^7r^p) = -cr{TT^7T^n) (35) 
8 4 

The quality of this relation is expected to improve with the dominance of the OPE graph. 
Of course, in the vicinity of threshold (p^ < SOOMeV/ c) one must take into account the 
isospin breaking (see Sec. 2. 4 above); in this case it would be more correct to consider the 
corresponding relations for quasi- amplitudes. 

Even a short glance at the Table 1 is enough to realize that the relations (^) strongly 
contradict the experimental data in the area of < 500MeV/c. As a consequence 
the OPE-dominance hypothesis is questionable and the straightforward application of 
the Chew-Low extrapolation procedure cannot result in reliable values of tttt parameters. 
This statement is now commonly agreed upon. 

II. The S'-wave final state. 

In accordance with the modern data on distributions [7-9,17,25] at < 380MeV/c 
each pair of final particles in reaction (|^) is expected to be in a relative S'-wave state. As 
it was shown in [^, in this case some of relations (^), namely 

a^TT^n'^p) = cr(7r^7r'^p) = -a(7r^7r^n) (36) 
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are expected to be valid. The comparison of ( p6D with the corresponding data from Table 1 
shows, however, the striking disagreement. Since both parts of the picture (distributions 
and total cross sections) contain the same dynamics we are forced to conclude that there 
is either an implicit contradiction between the data on distributions and those on total 
cross sections or that the assumption on S'-wave final states is too restrictive. 
It should be noted here that one of the relations ('B5l,B^): 



a(7r"7r°j9) = a{n+n^p) (37) 

seems to be supported by data. This relation implies (see Table 2) the smallness of the 
isotopic amplitude D in (^ in comparison with C. If this is true, all distributions in 
both channels must be approximately equal. From the purely theoretical point of view 
the smallness of D would imply a delicate cancellation of various contributions. That is 
why it would be extremely interesting to check this point at momenta > AOOMeV/ c 
where the existing data are still scarce. 

III. The phase-space like distributions. 

Before starting the discussion we would like to stress that our consideration here is 
only preliminary; as already stated in our introduction, the problems concerning the 
distributions will be discussed in detail in forthcoming publications. 

In accordance with the results of ref. 0, ||, |, ^ none of the distributions in 
four channels (except 7r~7r"'"n) up to < 320MeV/c show (within the error bars) any 
deviation from the purely 3-particle phase space kinematics. What does this fact tell us? 

To answer this question it is necessary to use eq. ( |A.8| ) along with the explicit expres- 



sions for the matrix elements Mij = 1, .., 4). A short look at eqs.(56) of ref. shows 
that almost each element of the matrix M depends, for example, on r. Thus to remove 
this ("necessary") r-dependence one has to select the suitable amplitude vectors X^. The 
direct calculation, however, shows that this is impossible: no physically reasonable choice 
of Xr can ensure the exact independence of IM^P of all four kinematical variables. This 
result becomes more clear from the examples considered in Appendix and the explicit 
form of the OPE-contribution. 

Thus, we are forced to conclude that the experimental observation of the "phase-space 
like" behavior of distributions would show that the precision of the data is insufficient to 
display the r - dependence originating from the OPE - graph. In fact it shows that the 
experiment cannot distinguish between spin 1/2 and spinless nucleons. Such precision, of 
course, does not provide the opportunity to extract reliable information on the details of 
the dynamics of the process under consideration. It can also result in large uncontrollable 
errors if the total cross section is obtained with the help of integration of distributions 
over the corresponding variables. 

So, even a very general consideration shows that the data set on total cross sections 
may contain insufficient information for the extraction of some parameters of the am- 
plitude. Moreover, some of the data points (especially, those obtained with the help of 



16 



the integration procedure) may be incorrect in addition to other - purely experimental - 
reasons. 



3.2 Results of the Fit 

In accordance with the general strategy of the data fitting as discussed in Sec. 2.6 we 
have performed 7 different fits, each one corresponding to the specific set of the amplitude 
parameters involved: B, n, A, B + n, B + A, vr + A, B + n + A. Some results of the 
best fits are given in Table 3. It should be noted here that for each kind of fit there were 
found several solutions with approximately the same values of x^] they differ from each 
other by the parameters values. Thus, in Table 3 we give the results corresponding to 
only one "typical" solution for each kind of fit. 

The results of the primary "filtration" of the data sample (exclusion of "doubtful data 
points", DP) show that the latter contains some DP's. Indeed, 7 points — namely, 6, 64, 
73, 75, 88, 91, 93 — appear in aU solutions as doubtful ones. The point no. 12 omitted 
from the very beginning must be also added to this list. These 8 points need special 
verification. Until the latter is not done, they must be excluded from the list of the data 
used for the subsequent analysis. Excluding them — the database is reduced by this 
procedure from 105 to 97 points — we can perform the second step of the data filtration 
to elucidate the status of 4 points {nn. 1, 4, 76, 90) which appear as doubtful in some 
solutions. The results of this — second — step are given in the Table 4. It is organized 
in the same manner as the Table 3; the points numbers placed in brackets indicate that 
the corresponding point gives the individual contribution xf close to the critical value 
ixl = 5) 

The consideration of Table 4 shows that two points - nn. 1,4 - which appeared in 
Table 3 as "almost doubtful" can be regarded as being really doubtful. For them we 
found no solution with xf < 4, {i = 1,4). 

Summarizing, the list of doubtful data points contains: 

Ti^Ti^n channel 1,4,6,12; 

7r°7r°n channel 64,73,75; (38) 

TT^TT^n channel 88, 91, 93, 

It should be noted that all DP's are concentrated in channels with the neutron in the 
final state; the two remaining channels (those with the final state proton) do not contain 
DP's. 

The database filtration procedure could be continued further, but we prefer to stop it 
here, because in the next paper we shall include the additional experimental information 



(the data of |^ on distributions in vr vr+n - channel) to get more reliable and well 



grounded results. It is necessary to note here that we do not regard the list (pSj) as the final 
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one. Unfortunately the "TRIUMF-OMICRON problem" [| in the 7r+7r+n-channel (which 
is clearly visible on fig.Se) is not as trivial as it seems at the first glance. The attentive 
reader could notice that the point no. 90 (TRIUMF) oftenly appears as the doubtful one 
along with the points nn. 88,91,93 (OMICRON) - at least in the best solutions. Our 
experience prompts that it is the hint on the impossibility to resolve the above mentioned 
problem with the help of purely statistical methods on the database in question: it is 
necessary to use additional independent information. As to our list of DP's in Tr^vr+n- 
and TT^TT^n-channels, we believe it is correct and the new information can only display 
some new DP's in addition to those quoted above. 

Our conclusion about the presence of DP's in the data sample in question contradicts 
the result derived by Burkhardt and Lowe (BL) in refs. [|l^, ^ where the authors have 
stated that the database on total cross sections is internally consistent up to piah = 
AOOMeV/c (apart from a momentum region 310 — 370MeV/c in the 7r"''7r"'"r;,-channel). It 
is interesting to elaborate on this aspect in more detail. 

In what follows we discuss only the latest results published in ref. |34|, where the new 
data |]2S| on total cross sections of the vr+Tr'^p-channel are taken into account. In accor- 
dance with the prescription of ref. all data points in Tr+vr+n-channel corresponding 
to the momentum interval piab = 310 — 370MeV/c (totally 7 points) must be excluded 
from the fitting procedure. Then, in their work Burkhardt and Lowe use the reduced 
data set: points nn. 4, 6, 7, 9, 11, 12, 16, 17, 41, 64, 70, 72, 74 are not included in 
their database. Thus, to compare the results, we have also excluded the points mentioned 
above and performed the fitting of the remaining database (49 points). Then, using the 
parameters corresponding to the best fit (or, if several, to one of them) we have computed 
the individual contributions xf ^or all points in the data set, including those previously 
omitted. The values of Xi in the BL fit, which are not listed in ref. have been ob- 



tained with the help of the interpolation of the numbers given in the Table 1 of ref. ||64 . 
The corresponding results concerning the DP's ( pSj ) are given in Table 5, where we use 
the parameters of one of the best solutions with the background amplitude B only (point 
no. 12, as noted above, is omitted from the very beginning). 

The comparison of numbers in Table 5 shows that the points in question look incon- 
sistent in both fits. This result confirms our list (|38|) of DP's. 



It must be noted also that in refs. |T3, points nn. 73, 75 have been omitted along 



with those quoted above, while near threshold points nn. 1, 3, 5 (and some others corre- 
sponding to higher momenta) were supplied with the beam momentum values consider- 
abley different from those quoted in the original papers, this latter difference originating 
from the additional corrections of purely experimental character which were taken into 
account in refs. |l^, |6^ and neglected in our analysis because we use only published re- 



1 The origin of the "TRIUMF-OMICRON problem" has been recently pointed out by Ortner et 
al. For the OMICRON collaboration the total cross sections have been obtained from angular 

distribution data in a restricted angular range assuming a smooth continuation over the full phase space. 
Angular distributions for the n^p — > TT^7T~n channel um do not support such a simple extrapolation. 
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suits. (We are grateful to Prof. J.Lowe for a detailed discussion of this question). In 
addition, the systematic errors in the vr^Tr^n-channel were handled in a different way in 



refs. [0, The differences pointed out above between the two data samples give rise 
to the different corresponding conclusions. 

It should be mentioned that if the corrected values of effective beam momentum (in 
accordance with BL data base) are used in our analysis, only the point no. 1 (corre- 
sponding to = 302.5 MeV/c in the BL sample) looks considerably better and must be 
deleted from the list (|38D. All other results remain practically unchanged. 

From our findings above a detailed discussion of the parameter dependence of the 
model would be premature. This can be best illustrated by the mutual comparison of 
numbers from Table 6 where we give the typical examples of solutions obtained in a course 
of the pruned database fitting with different groups of parameters. The considerations of 
this table shows that every statement on the values of TTTT-parameters based on the analysis 
of the data on total cross sections only would be strongly model dependent because the 
present knowledge of the background processes (and, hence, the values of Ai, . . . , An) is 
too pure, the latter processes playing a very important role in the nN tttcN reaction 
mechanism. 



4 Conclusion 

Before the discussion of our results we would like to start with some notes concerning the 
adequate treatment of the experimental data on vrA^ — >■ nirN reaction in the momentum 
region in question (p^ < 500MeV/c). 

Even a short glance at Figs. 3a-3e is enough to recognize that the modern database 
contains some pairs of mutually inconsistent points. Thus, the preliminary selection (or, 
the same, filtration) of data is unavoidable if one intends to extract with the reasonable 
accuracy such quantities as vrvr-scattering parameters. The only question is how to do 
this selection. 

Of course, one can simply exclude the mutually contradictory pairs of points from the 
database as it was done in refs. [^, |6^ . Unfortunately, such a recipe does not guarantee 



the consistency of the remaining collection of data. Indeed, after removing "doubtful 
pairs" one still has 5 independent subsets of data points (the total cross sections of 5 
channels). These are related by isospin and Bose symmetries, and it is essential that 
these symmetries are incorporated in the analysis. The channels are also related by 
crossing symmetry (which was ignored in the analysis of refs. [ffR 0)' 



Thus, it is evident that any correct data selection must take into account all symmetry 
requirements. If this condition is fulfilled one can use a wide range of equivalent approx- 
imants for the "data filtration" . Their specific choice (polynomials, spherical harmonics, 
etc.) can have influence only on the number of free parameters necessary to achieve a 
good description of the data. From this point of view our choice [B + U + A) for the 
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vrA^ — s> vrvrA^ transition amplitude seems to be the most economic one. In addition to 
the symmetry requirements it takes account of rather general features of the process in 
question: the existence of a smooth background B and the vr and A poles close to the 
kinematical regime investigated. 

Of course, we could use polynomials of higher degree for the background contribution 
(this would introduce ~ 30 new parameters for the second degree polynomials). Our 
results show, however, that this is not necessary, because a reasonably good description 
of data is achieved with linear approximants (in our next paper it will be shown that this 
is true also for the data on distributions). 

Summarizing we conclude that the procedure used here for the primary filtration of 
data is general enough: in the case when some data points cannot be described with the 
reasonable accuracy, they need further experimental verification. Until such a verification 
is not done these points must be excluded from the database. 

With this in mind this conclusion we can make now the following statements. 

• The modern collection of the experimental data on total cross sections of 5 channels 
of nN — > TiiiN reaction at < 500MeV/ c contains some doubtful points (as listed 
in eq. (|38|) ) which need special verification. 

• The experimental data on two channels with a final state proton {tt^tt^p, tt^tt^p) 
are in statistical sense reliable: they do not contain DP's. Nevertheless, it would be 
extremely useful to perform new accurate measurements of the 'K^Ti^p channel to 
increase its statistical significance. It would be of special interest to perform such 
measurements at relatively large values of the incident momentum: = 400 — 
500MeV^/c. The corresponding results would give the opportunity to check the 
approximate equality of cross sections of these two channels. 

• The pruned database (i.e. that obtained after the exclusion of DP's) can be de- 
scribed in terms of our ansatz with a reasonable value of (< 1-2) (see Fig. 3). 
Unfortunately, there exist at least several quantitatively equivalent solutions corre- 
sponding to rather different values of parameters (the correlations among them are 
extremely strong). Thus, the experimental data in question do not allow the extrac- 
tion of a unique set of parameters in the ttN ttttN. Only qualitative conclusions 
concerning the influence of the various reaction mechanisms are possible from this 
limited set of data. 

• The most important contribution comes from the background term B. This term 
alone admits solution with = 1.4 The addition of the 11 term improves the de- 
scription significantly: x^{B + 11) = 1.2. This fact shows that — in accordance 
with common expectations — the OPE term plays an essential role, its contri- 
bution is definitely visible even on the database in question and, therefore, the 
corresponding parameters might be extracted if the database is enlarged by the in- 
clusion of the accurate enough data on distributions. The role of 11 term, however. 
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is by no means dominant (x^(n) = 2.3), so the well known Chew-Low extrapola- 
tion procedure seems questionable. This finding coincides well with the results of 
refs. [||,|5|, H, g. 

• In contrast to the 11 term, the A contribution is rather insensitive on the database 
in question. The only observation which is of some interest can be formulated as 
follows: the most significant A contibution corresponds to the parameters Qq and 

appearing in graphs of Fig. 2d. No definite statements concerning the other 
four parameters associated with the mechanism of the isobar exchanges (graphs of 
Figs. 2a-2d) can be done. In particular, the values of g-j^NA and g^AA cannot be 
fixed from the data on total cross sections. 

To study the A mechanism one needs the more extensive database which must 
contain the data on distributions at several beam momenta. 

• The informative capacity of the modern database on total cross sections is absolutely 
insufficient for the model independent extraction of the parameters describing the 
low energy vrvr interaction. The widely used statement that these parameters are 
most sensitive on and, thus, best extracted from the data on total cross sections near 
threshold, which is based solely on the Olsson-Turner model, is misleading. For the 
extraction of reliable information on vrvr interaction it seems much more reasonable 
to use the accurate data on distributions at several values of incident momentum in 
the interval 300 - 500 MeV/c (see ref. 0). 

Summarizing the present status of our analysis we conclude that information from the 
total irN — s> vrvrA^ cross sections only rather quantitatively fixes the basic features of the 
dynamics. Evidently more sensitive data, such as angular distributions, are crucial for a 
quantitative understanding. The analysis of the enlarged database containing data both 
on total cross sections of 5 channels and on distributions in the 7r~^7T~n channel is presently 
still troubled with serious numerical problems. It will be discussed in a forthcoming paper. 
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A Appendix: "Quasi" nN elastic scattering 

As an illustrative example on the subtle influence of different dynamical elements on the 
extrapolation and the angular distributions in different isospin channels, let us consider 
the familiar process: "quasi" nN elastic scattering: 

Ti'^i^k) + N:{p) ^ Tx\k') + N>}{p') (A.l) 

The term "quasi" used above means that — to simplify the reaction kinematics — we 
consider the hypothetical case of equal pion and nucleon masses k = p = k' = p' = m 
In ( |A.1| ) a,b = 1,2,3 and e, / = 1,2 are the isotopic indices, /i, z/ = 1,2 the nucleon 
polarizations. 

The amplitude of the process ( |A.1|) can be represented in the following form (for a 
comprehensive review see ref. ||61|| ): 



TZfe = SbaSfe T^'^'-^^ + Z6tac{rc) fe T'^^^"^ (A.2) 

Each of the isotopic amplitudes T^^'^^'^^ can be written as follows: 

2-mK±) = u^{p'){mA^"'^^'> + QB^"'^^'>}u''{p) (A.3) 

here Q = {k + k')/2, A"^ and 5^ being the functions of scalar kinematical variables: 
t and u = (s — m)/4. The coefficient m at A"^ in ( |A.3| ) is introduced for dimensional 
reasons. The crossing and Bose symmetries imply: 

A^{u,t) = ±A^{-u,t) 

B^{u,t) = TB^{-u,t) (A.4) 

Formulae (|A.2| - |A.^) determine the so-called "canonical" form of the amplitude Tj^afe 
of the process ( [A.ll ). 

Let us consider now three experimentally accessible channels of the reaction ([A. 1|) and 
write down the corresponding amplitudes (T = A,B): 

1. n+p^n+p Ti = T+-^T- 

2. 7i-p^7i-p T2 = T+ + ^T- (A.5) 

3. n^p — s> n^n T3 = — -^T~ 

Evidently, the amplitude of the r-th channel can be written in the form of eq. ( [A.3| ) 
with Ar,Br being constructed in accordance with eqs. (|A.5|) . 
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Suppose, further, that we analyze the data on three channels ( |A.5| ) with unpolarized 



nucleons. In this case the amplitude modulo squared (summed over the final and averaged 
over the initial nucleon polarizations) has the form: 

|T^|2 = ^Sp{{p + m){mA^ + QB^){p + m){mA + QB)} (A.6) 

or, equivalently: 

\Tr\'= E Vk^Gk^V^ (A.7) 

k,p=l,2 

here, = {A'^,B''') is the "amplitude vector" of the r-th channel and the real and 
symmetric matrix G has the following elements: 

Gil = m\4m^ - t), Gu = G21 = 4m^u, G22 = + ^t{4m^ - t) (A.8) 
Thus, the differential cross sections for the r-th channel may be written as follows: 

dar = N ■ V^GV ■ dt (A.9) 

where N contains all "trivial" kinematical and normalization factors. 

Let us suppose that it is possible to make use of the "smoothness hypothesis": the 
amplitudes and B^ in the low energy region area are smooth functions of their variables 
z/, t: 

yl+(z/,t) = ao + ai ■ t A~{i',t) = ^av 
B-{u,t) = 2{(3q + (3i-t) B+{u,t)=Ahu (A. 10) 



The expressions ( [A. 101 ) are the most general linear forms consistent with the require- 



ments ( |A.4|) . These forms are determined by 6 parameters (for the sake of simplicity we 
consider them to be real) which describe the cross section of each channel ( |A.5D . 

Let us consider now 4 instructive examples I,..., IV for the sensitive interplay of the six 
parameters. In each case we write down explicitly: the values of the parameters (a), the 
exact expressions for |Trp in terms of s,t (b) and the "extrapolation laws" for the total 
cross sections 0"^. ~ / iT^prft, i.e. the leading order terms of the low energy expansion of 

in power series of CMS momentum fc^ (c), together with short comments (d). 
Example I 

a. ao = cti = a = & = /?i = 0, Po = d 

b. iTip = IT2I2 = llTgp = d^{{s - 2m2)2 + (s - m^)t} 

c. ai ~ a2 ~ |cr3 ~ A;°(l -h 0(F)) 

d. The angular distributions are relatively simple (linear in t) in all channels, moreover, 
they are identical due to the presence of the dynamical symmetry (i.e. conditioned solely 
by the particular values of the parameters). The form of the extrapolation law in each 
channel is "natural" (~ k^) 
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Example II 

a. ao = ai = b = Po = Pi — 0, a — d/Am^ 

b. \T,\' = \T,\' = l\Ts\' = 

c^Vl6m2{16m2(s - 2m'^)^ - 4(s - 2m'^){s - 6m'^)t - 4(s - 3m'^)t^ - t^} 

C. (7i~(72~icT3~fc0(l + O(A;2)) 

d. The same situation as in I with the only difference, that the angular distributions 
arc more complicated. 
Example III 

a. ai = a = b = Pi = 0, ao = Po = d 

b. \Ti\^ = d^{{s - Am^f + {s - 4:m^)t} 
\T2\^ ^ d''{s^ + St} 

iTsp = 2d^{{s - 2m2)2 + (s - m^)t} 

c. (Ti ~ A;^ < a2 ~ 2(T3 ~ 2m4(l + 0{k^)) 

d. Simple angular distributions (linear in t) in all channels. The "unnatural" form of 
the extrapolation law (~ k'^) in the first channel, opposite to the natural forms in the 
third and second channels. In addition here is a dynamical symmetry relation between 
the total cross sections (T2 and (T3. 

Example IV 

a. ai = a = b = Pi = d/lGni^, ao = Po = —d/4 

b. |Ti|2 = (i7(16m2)2 . 4s2{(s - Am^f + (s - Am^)t} 

iTsp = dV(16m2)2 . 4s{(s - Am'")'' + {s - 4m2)(3s - 4:m^)t 
+ (3s-8m2)t2 + t3| 

iTsp = ci7(16m2)2 . 2{-4m2s(s - Am'^)t + s{s - 9>m^)t^ - mH^} 

c. (Ji ~ 3/2(T2 — 60-3 ^ 

d. The angular distributions in three channels are quite different, the term ~ t° being 
absent in the third channel. Unnatural forms (~ k^) of the extrapolation law show up in 
all channels. 

Table Captions 

Table 1: 

List of all data for the total cross section from threshold up to k^i'^ = 500 MeV/c for 
various nnN final states. 

a) The 7r''"7r~n- channel 

b) The TT^TT^p- channel 

c) The 7r°7r°n - channel 

d) The 7r+7r°p- channel 

e) The tt+tt+ti - channel 

Table 2: 

Basic channels of the ttN —>■ ttttN reaction. The isospin amplitudes A,B,C,D are defined 
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in eq. (|). 



Table 3: 

The results of the first step of the data filtration. The different columnes denote: 
Column 1: the kind of the fit. 

Column 2: the total value of corresponding to one of the best solutions. 
Column 3: the value of (from eq. (pl|)) 

Column 4: the list of the "doubtful points" (DP's); (i.e. points Pq with > 5) numbers 
are given in accordance with Table 1. 

Column 5: summed contribution x^{DP) of all DP's listed in column 4 to the total value 
of x^- 



Table 4: 

The results of the second step of the data filtration. The notation is as in Table 3, the 
numbers in bracket refer to points with a close to 5. 



Table 5: 

Comparison of the for selected data points from different ttttN channels in the analysis 
from Burkhardt and Lowe ||T8| , |64| with the results of this analysis. 



Table 6: 

Examples of solutions obtained in a course of the fitting with different groups of param- 
eters. 
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Table 1: List of all data for the total cross section from threshold up to k'f'^ = 500 MeV/c 
for various nnN final states, 
a) The vr+Tr^n- channel 



channel 


^ in MeV/c 


atot in mb 


error 


Mi 


error 


reference 
















1 


295 


0.0051 


0.0012 


862 


101 




i 




2 


313 


0.0138 


0.0015 


685 


37 


[ 






3 


315 


0.020 


0.0031 


781 


60 




i 




4 


321 


0.015 


0.003 


584 


58 


[ 


42] 


5 


334 


0.051 


0.012 


837 


99 




i 




6 


334 


0.027 


0.005 


609 


56 


[ 


42 




7 


336 


0.031 


0.02 


624 


208 




8 


342 


0.0603 


0.0032 


801 


21 


[52] 


9 


346 


0.053 


0.013 


712 


87 


[ 






10 


354 


0.118 


0.020 


959 


81 




i 




11 


360 


0.125 


0.028 


922 


103 




42 




12 


361 


0.060 


0.015 


638 


80 


[54 




13 


369 


0.166 


0.0061 


965 


18 


[52] 


14 


374 


0.211 


0.036 


1039 


89 




i 




15 


374 


0.14 


0.10 


843 


301 


[55 




16 


378 


0.16 


0.06 


869 


163 


i 


42 [ 


17 


392 


0.4 


0.2 


1229 


307 


m 


18 


394 


0.327 


0.041 


1095 


69 




i 




19 


395 


0.374 


0.0146 


1161 


23 


[52] 


20 


404 


0.38 


0.09 


1098 


130 


[ 


42 




21 


408 


0.546 


0.0306 


1279 


36 


[52 




22 


413 


0.477 


0.056 


1163 


68 




i 




23 


415 


0.57 


0.06 


1253 


66 


[ 


4f:] 


24 


428 


0.71 


0.17 


1333 


88 


[55 




25 


432 


0.785 


0.104 


1249 


150 




i 




26 


449 


1.160 


0.052 


1483 


33 




52 




27 


451 


1.052 


0.125 


1409 


84 




i 




28 


456 


1.0 


0.2 


1338 


134 


[ 


44 




29 


458 


1.39 


0.05 


1562 


28 


[ 


43 




30 


477 


1.880 


0.077 


1677 


34 


[52[ 


31 


485 


1.93 


0.16 


1645 


68 




IC 




32 


485 


2.4 


0.26 


1834 


99 




41 




33 


491 


1.93 


0.37 


1608 


154 


[55 




34 


495 


2.6 


0.26 


1840 


92 




41 




35 


499 


2.40 


0.16 


1740 


58 




44 





29 



Table 1: b) The tt vr^p- channel 



channel 


j^iab MeV/c 


(Xtot in mb 


error 


Mi 


error 


reference 


7T~7T°p 














36 


295 


0.00075 


0.0004 


216 


58 







37 


315 


0.0022 


0.0007 


210 


33 


[0 




38 


334 


0.0085 


0.0016 


301 


28 







39 


354 


0.020 


0.005 


356 


44 


[0 




40 


375 


0.027 


0.006 


342 


38 







41 


391 


0.08 


0.13 


517 


420 


H 


1 


42 


394 


0.050 


0.013 


400 


52 


m 




43 


406 


0.2 


0.1 


738 


184 


wn 


44 


413 


0.073 


0.015 


428 


44 


m 




45 


415 


0.11 


0.05 


520 


118 


|47| 


46 


415 


0.09 


0.01 


469 


26 


|46 




47 


428 


0.14 


0.07 


526 


142 


|41| 


48 


432 


0.119 


0.020 


493 


41 


[0 




49 


450 


0.157 


0.037 


520 


61 


m 




50 


456 


0.17 


0.05 


527 


78 


P41 


51 


458 


0.17 


0.01 


523 


15 


|43| 


52 


495 


0.31 


0.07 


613 


69 


|41| 


53 


499 


0.32 


0.05 


613 


48 


|44| 
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Table 1: c) The - channel 



channel 



f^iab MeV/c 



atot in rnh 



error 



M, 



error 



54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
66 
67 
68 
69 
70 
71 
72 
73 
74 
75 
76 



272 
276 
280 
284 
286 
287 
291 
293 
298 
305 
310 
314 
323 
331 
339 
349 
353 
359 
385 
390 
391 
400 
494 



0.00038 

0.00059 

0.00118 

0.00206 

0.00231 

0.00333 

0.00381 

0.0081 

0.0085 

0.0171 

0.032 

0.0219 

0.0303 

0.0598 

0.0752 

0.0981 

0.13 

0.118 

0.32 

0.388 

0.27 

0.479 

1.3 



0.00010 

0.00014 

0.00022 

0.00035 

0.00065 

0.00064 

0.00081 

0.0013 

0.0010 

0.0019 

0.005 

0.0020 

0.0030 

0.0064 

0.0073 

0.0093 

0.02 

0.011 

0.04 

0.046 

0.07 

0.049 

0.1 



535 

466 

465 

477 

461 

523 

471 

648 

563 

661 

810 

615 

620 

772 

770 

785 

870 

781 

1040 

1108 

913 

1151 

1236 



70 

55 

43 

40 

65 

50 

50 

52 

33 

37 

63 

28 

31 

41 

37 

37 

67 

36 

65 

66 

118 

59 

48 
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Table 1: d) The tt+tt"]?- channel 



channel 


yab MeV/c 


cTtot in mb 


error 


M, 


error 


reference 
















77 


298 


0.001 


0.0017 


223 


190 






1 


78 


309 


0.0027 


0.0012 


267 


59 


i 


m 


79 


331 


0.0068 


0.0018 


278 


37 


I 


251 


80 


342 


0.018 


0.012 


387 


129 


|48 




81 


352 


0.0146 


0.0026 


309 


28 


i 


25 


1 


82 


374 


0.026 


0.0027 


336 


17 




83 


392 


0.048 


0.034 


399 


141 


|48 




84 


415 


0.12 


0.05 


542 


113 


|47| 


85 


418 


0.189 


0.036 


669 


64 




i 
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Table 1: e) The vr+vr+n- channel 



channel 


j^iab MeV/c 


atot in mb 


error 


M, 


error 


reference 


TT+TT+n 














86 


288 


0.00011 


0.00004 


236 


32 


m 


87 


292 


0.00028 


0.00005 


249 


22 


1 


m 


88 


298 


0.0018 


0.0004 


444 


49 








89 


299 


0.00060 


0.00010 


243 


20 


1 


16] 


90 


310 


0.00146 


0.00022 


246 


19 


1 


16] 


91 


317 


0.0080 


0.0018 


467 


53 




P 




92 


335 


0.0094 


0.0023 


338 


41 


[50 




93 


338 


0.0217 


0.0045 


513 


53 








94 


342 


0.030 


0.018 


565 


170 


1 


48 




95 


358 


0.0274 


0.0052 


440 


42 








96 


364 


0.0228 


0.0048 


378 


40 




5C 




97 


378 


0.039 


0.007 


432 


39 




i 




98 


390 


0.026 


0.055 


317 


336 


1 


48 




99 


398 


0.0451 


0.0103 


395 


45 




i 




100 


418 


0.065 


0.0135 


415 


43 




i 




101 


430 


0.0537 


0.010 


353 


33 


[5C] 


102 


439 


0.074 


0.0153 


394 


41 




i 




103 


459 


0.083 


0.0178 


380 


41 




i 




104 


477 


0.12 


0.01 


424 


18 


[38 




105 


480 


0.094 


0.0200 


370 


39 




i 
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Table 2: Basic channels of the vrA^ rnvN reaction. The isospin amplitudes A, B, C, D 
are defined in eq. (||). 



channel 




amplitude 


st at. fact or / 


1T~p —>■ TC'TT^n 


1 


Xi = fiA + C) 


1 




2 


X2 = \C-D 
X, = ^A 


1 


Ti^p —>■ TT^n^p 


3 


1/2 


Ti+p 7r+7r°p 


4 


X^ = \C + D 


1 




5 


= ^(B + C) 


1/2 
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Table 3: The results of the first step of the data filtration. The different columnes denote: 
Column 1: the kind of the fit. 

Column 2: the total value of corresponding to one of the best solutions. 
Column 3: the value of ^ (from eq. (^)) 

Column 4: the list of the "doubtful points" (DP's); (i.e. points Pq with Xg ^ 5); numbers 
are given in accordance with Table 1 

Column 5: summed contribution x^i^P) of all DP's listed in column 4 to the total value 
of x^- 



Fit 






List of DP's 


x'iDP) 


B 


188.6 


2.03 


4,6,64,73,75,76,88,91,93,95 


87.6 


TT 


284.7 


2.85 


many 




A 


1396.0 


14.2 


many 




B + n 


149.0 


1.66 


1,6,64,73,75,88,90,91,93 


64.6 


B + A 


180.8 


2.08 


4,6,64,73,75,76,88,91,93 


78.1 


71 + A 


204.4 


2.17 


4,6,64,67,72,73,75,76,85 ,88,90,91,93,95 


109.7 


B + 71 + A 


148.2 


1.79 


1,6,64,73,75,88,90,91,93 


64.1 
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Table 4: The results of the second step of the data filtration. The notation is as in Table 
3, the numbers in bracket refer to points with a close to 5. 



Fit 




x' 


List of DP's 


X\DP) 


B 


123.4 


1.43 


[1]4,67,72,76,95 


30.6 


TT 


216.3 


2.33 


many 




A 


1259.0 


13.8 


many 




B + n 


98.8 


1.20 


1,[4],90 


11.8 


B + A 


119.1 


1.49 


[1],4,67,72 


17.3 


71 + A 


139.5 


1.60 


4,19,21,67,72,85,95,97 


52.6 


B + 71 + A 


98.2 


1.29 


1,[4],85 


12.0 



Figure Captions 

Figure 1: 

Schematical decomposition of the nN ttttN reaction amplitude into the background 
part {B), the A- isobar piece (A) and the OPE-term (11). 



Figure 2: 

Different contributions of the A-isobar. (a), (b), (c): A-excitation in the first, second and 
both intermediate states, respectively; (d): nonlinear ttttNA coupling. 



Figure 3: 

Comparison of the model independent fit to the quasi amplitude Mi (eq. (pOD) in 5 
vrvrA^ channels with the data as a function of p-^ (GeV/c). Note that "doubtful points", 
as discussed in Sect. 3.4, are not included in the data. The experimental data and 
corresponding references are from Table 1. Data points from Kernel et al. P, 0, ^ are 
emphasized by asterisks. 

a) The 7r"*"7r~n- channel. 

b) The TT^vr^p- channel. 

c) The ir^TT^n- channel. 

d) The 7r"'"7r°p- channel. 
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Table 5: Comparison of the for selected data points from different vrvrA^ channels in 
the analysis from Burkhardt and Lowe |18, 64] with the results of this analysis. 



no. point 


Individual contribution xf 


(channel) 


BL 


Our results 


1 {iT^TT^n) 


5.8 


5.8 


4 (TT^TT+n) 


4.2 


3.9 


6 (7r~7r"'"n) 


9.6 


9.1 


64 (tt^tt^) 


8.1 


8.4 


73 (7r°7rOn) 


4.4 


3.7 


75 {TT^ir^n) 


5.3 


3.5 


88 (vr+TT+n) 


10.0 


9.0 


91 (vr+TT+n) 


8.9 


7.4 


93 (vr+TT+n) 


10.4 


7.7 



e) The Tr+Tr+n- channel. Compared are the fits with the full data set and with a partial 
data set (lower and upper line, respectively, 4 points out of 20 are suppressed). 
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Table 6: Examples of solutions obtained in a course of the fitting with different groups of 
parameters 



B{x^ = 1.4) 




2.3) 


S + 7r(x' = 2.3) 


Para- 


value 


error 


correlation 


value 


error 


correlation 


value 


error 


correlation 


meter 


10^ 


103 


factor 10^ 






factor 10^ 


103 


10^ 


factor 10"^ 




0.39 


1.12 


4.48 


0.0 


fixed 




-2.70 


0.18 


0.20 


A2 


-5.00 


10.9 


6.79 


0.0 


fixed 




-4.46 


4.58 


2.37 


A, 


-10.4 


36.7 


27.2 


0.0 


fixed 




-12.9 


7.56 


2.87 




-0.77 


0.35 


2.11 


0.0 


fixed 




-2.10 


0.34 


10.6 


A, 


3.77 


11.1 


21.2 


0.0 


fixed 




-4.43 


6.59 


21.8 


A, 


0.76 


6.09 


25.8 


0.0 


fixed 




5.63 


0.34 


0.30 


A, 


-6.11 


7.67 


4.43 


0.0 


fixed 




-3.16 


11.6 


48.9 




0.85 


1.86 


3.66 


0.0 


fixed 




1.84 


1.81 


17.9 


^9 


6.34 


6.19 


1.09 


0.0 


fixed 




7.88 


1.13 


0.39 


AiO 


0.41 


0.55 


0.26 


0.0 


fixed 




-0.15 


0.11 


0.06 


All 


8.86 


6.33 


0.44 


0.0 


fixed 




0.37 


2.14 


0.40 


Go 


0.00 


fixed 




9.7 


3.53 


4.7 


-0.18 


0.03 


0.57 


Gi 


0.00 


fixed 




-672 


44.4 


3.1 


-2.74 


0.17 


0.12 


G2 


0.00 


fixed 




1.5 • 10^ 


427 


4.7 


20.4 


12.2 


3.36 


Gs 


0.00 


fixed 




2.9 • 10^ 


405 


3.1 


18.7 


11.3 


2.56 
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Figure 2: Different contributions of the A-isobar. (a), (b), (c): A-cxcitation in the 
first, second and both intermediate states, respectively; (d): direct nnNA coupling (the 
crossing graphs are implied to be taken into account also). 
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Figure 3: Comparison of the model independent fit to the quasi amphtude Mj (eq. (|30D ) in 
5 uttN channels with the data as a function of {GeV/c). Note that "doubtful points", 
as discussed in Sect. 3.4, are not included in the data. The experimental data and 
corresponding references are from Table 1. Data points from Kernel et al. P, 0, ||, ^ are 
emphasized by asterisks. 

a) The vr+vr^n- channel. 
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Figure 3: b) The tt ttV" channel. 
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Figure 3: c) The Tr'^Tr^n- channel. 
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Figure 3: d) The tt'^ttV" channel. 




Pn iGeV/c) 
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Figure 3: e) The 7r+7r+n - channel. 

Compared are the fits with the fuh data set and with a partial data set (lower and upper 
line, respectively, 4 points out of 20 are suppressed). 
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